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Abstract 

We extend previous analyses of the supersymmetric loop correction to the neutral 
Higgs couplings to include the coupling XjXkHf. The analysis completes the pre- 
vious analyses where similar corrections were computed for the frHf, bbH®, ccH® 
and for itHf couplings within the minimal supersymmetric standard model. The 
effective one loop Lagrangian is then applied to the computation of the neutral 
Higgs decays. The sizes of the supersymmetric loop corrections of the neutral 
Higgs decay widths into xtxj {i — 1)2; j = 1,2) are investigated and the super- 
symmetric loop correction is found to be in the range of 7 ~ 15% in significant 
regions of the parameter space. By including the loop corrections of the other 
decay channels bb, it, fr, cc, and x°iX°j (* = 1 ~~ 4; j = 1 — 4), the corrections to 
branching ratios for — > xtxj can reach as high as 40%. The effects of CP 
phases on the branching ratio are also investigated. 



Permanent address. 



1 INTRODUCTION 



The neutral Higgs couplings to different fields are of great current interest as they 
enter in a variety of phenomena which are testable in low energy processes [1]. It is 
known that supersymmetric corrections can affect the neutral Higgs boson decays 
into bb, rf and cc. The decay properties of the lightest Higgs boson in MSSM would 
be different from those of the Standard Model Higgs boson when these corrections 
are taken into consideration. Specifically the ratio of the branching ratios to bb 
and rf of the Higgs boson is an important piece of evidence that might distinguish 
between the lightest MSSM Higgs boson and the Standard Model one at colliders. 
In MSSM there are also other modes for neutral Higgs decays that do not exist in 
Standard Model such as charginos and neutralinos. 

In this paper we compute the one loop corrected effective Lagrangian for the 
neutral Higgs and chargino couplings. We then analyze the effects of the loop 
corrections to the neutral Higgs decays Hf — > XjXk- m the analysis we also 
include the effect of CP phases arising from the soft SUSY breaking parameters. 
It is well known that large CP phases can be made compatible [2, 3, 4] with 
experimental constraints on the electric dipole moments (edms) of the electron [5], 
of the neutron [6], and of the Hg 199 [7]. Further, if the phases are large they could 
affect the Higgs sector physics. It is well known that one loop contributions to 
the Higgs masses from the stop, sbottom, the chargino and neutralino sectors can 
lift the lightest Higgs mass above M z . The inclusion of the CP violating phases 
brings mixings between the CP even and the CP odd Higgs [8, 9, 10, 22, 23, 24]. 
The CP violating phases modifies the physics of dark matter [11], and of other 
phenomena [12]. (For a review see Ref.[13].) 

The current analysis of A£#o x + x - and neutral Higgs decay into charginos is 
based on the effective Lagrangian method where the couplings of the electroweak 
eigen states H\ and if| with charginos are radiatively corrected using the zero 
external momentum approximation. The same technique has been used in calcu- 
lating the effective Lagrangian and decays of Hf into quarks and leptons [1, 15, 16]. 
It has been used also in the analysis of the effective Lagrangian of charged Higgs 
with quarks [1, 17] and their decays into tb and u t t [18] and into chargino + neu- 
tralino [19]. The neutral Higgs decays into charginos have been investigated before 
in the CP conserving case [20, 21]. In these analyses, the wave function renormal- 
ization and the counter terms for the mass matrix elements are calculated beside 
the vertex corrections of the mass eigen states h°, H° and A with charginos. In 
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the effective Lagrangian technique with zero external momentum approximation, 
the radiative corrections of the processes considered here originate only from the 
vertex contributions. Thus our analysis of the neutral Higgs decays into charginos 
is a partial one. However, as mentioned before the above analyses were carried out 
in the CP conserving scenario. As far as we know, the analysis for the neutral Higgs 
decays into charginos, with one loop corrections, in the CP violating case where 
the neutral Higgs sector is modified in couplings, spectrum and mixings, does not 
exist. We evaluate the radiative corrections to the Higgs boson masses and mixngs 
by using the effective potential approximation. We include the corrections from 
the top and bottom quarks and squarks [22], from the chargino, the W and the 
charged Higgs sector [23] and from the neutralino, Z boson, and the neutral Higgs 
bosons [24]. It is important to notice that the corrections to the Higgs effective 
potential from the different sectors mentioned above are all one-loop corrections. 
The corrections of the interaction AjC h o x + x - to be considered in this work are all 
one-loop level ones. So the analysis presented here is a consistent one loop study. 

The outline of the rest of the paper is as follows: In Sec. 2 we compute the 
effective Lagrangian for the XjXk-Hi interaction. In Sec. 3 we give an analysis 
of the decay widths of the neutral Higgs bosons into charginos using the effective 
Lagrangian. In Sec. 4 we give a numerical analysis of the size of the loop effects 
on the partial decay width and on the branching ratios. Conclusions are given in 
Sec. 5. 

2 LOOP CORRECTIONS TO NEUTRAL HIGGS 
COUPLINGS 

The tree-level Lagrangian for XjXk~H° interaction is 

C = ^xfPnxiHl + ^^PnxtHl + H.c, (1) 

where Hi and H\ are the neutral states of the two Higgs isodoublets in the minimal 
super symmetric standard model (MSSM), i.e., 




and the couplings 0^ and ipjk are given by 

4>jk = -gUwVj!, ip jk = -gU kl V j2 (3) 
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where U and V diagonalize the chargino mass matrix so that 



U*M X + V" 1 = diag(m x +, m x +) (4) 

The loop corrections produce shifts in the couplings of Eq. (1) and the effective 
Lagrangian with loop corrected couplings is given by 

(i/> jk + 5ij jk )4P R xtH 2 2 + A^ k xfP LX tHl + H.c. (5) 

In this work we calculate the loop correction to the XjXk-^i usm g the zero external 
momentum approximation. 

2.1 Loop analysis of 5(f)jk and Aipj^ 

Contributions to 5(f) j k and Aif)j k arise from the thirteen loop diagram of Fig. 1. 
We note that the contribution from diagrams which have H + W + H°&nd H°Z°H° 
vertices do not contribute in the effective Lagrangian with zero external momentum 
approximation since these vertices are proportional to the external momentum. We 
discuss now in detail the contribution of each of these diagrams in Fig. 1. We 
begin with the loop diagram of Fig. li(a) which contributes to 5<f>jk and Aipj k . 
We calculate the corrections of the amplitude from Fig. li(a) 

SM = idcpjkUjPRVk + iAtp jk UjP L v k (6) 

The idea is to extract, from the amplitude correction, the expressions for 5(fij k and 
Aipj k from those parts that are proportional to UjPnv k and UjPiv k respectively. 
For this purpose we need bbH\ interaction which is given by 



Ash? = HubWHl + H.c. (7) 



where Hu is given by 



H a = - {{-\ + \ sin 2 e w )D* m D m - \ sin 2 6 w D* b2l D b2l ) cos/3 

V 2 cos i) w Z 6 6 



{Dl u D m + D* m D b2l ) - D* b2l D bll (8) 



y/2mw cos (3 1 1 y/2mw cos (3 



The matrix elements D q are defined as 



D+MlD q = diag{m%,m%) (9) 
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We need also the t\ + b interaction which is given by 

Ax + B = -9xt[(U* kl D* bli - K b U* k2 DljP L 

-K t V k2 D* bu P R ]tb* + H.c (10) 

where K ti b are given by 



\/2m w sin (5 
V 2m w cos p 

For external momenta s, q and q — s the amplitude correction from loop li(a) is 
given by 

SM = -g 2 H u u(q - s)[C Ljl P L + C Rjl P R ] 

j f)+^][ci k p R +c Rk p L Hs) 

X ((s + £) 2 -m 2 + ie)(£ 2 - ml + ie)((£ + qf - m| + ie) (12) 
where Ci j{ and C Rjl are given by 

c Ljl = u^di - K b u; 2 D* b2i 

C Rjl = -KtV j2 D* hu (13) 

The part in the numerator 

[C Ljl P L + C Rjl P R }y+f)+m t } 

{C* Lm Pr + C* Rki P L ) (14) 

could be written as 

+m t [C Rjl C* Lki P R + C Ljl C Rki P L \ (15) 

by using the facts that ^P L = P R ^, P L P R = 0, Pf = P L and P R = P R . The 
first term in Eq. (15) does not contribute to S(pjk or Aipj k since it does not have 
the same Lorentz structure. The second term of Eq. (15) contributes the part of 
m tCR jt Cl k . to 5(pjk and ratC^C^ . to Aipj k . Thus the loop corrections 5<pjk and 
Aipjk read 

iS(f) jk = -g 2 H a m t C R]l C* Lk J 
iAip jk = -g 2 H il m t C L , jl C Rk J 

(16) 
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where 

r _dH 1 

J (2tt) 4 ((s + t) 2 -m 2 t + ie)(£ 2 - m| + ie)((£ + q) 2 - m| + ie) { } 

Now for zero external momentum approximation we set s = q = 0, and the integral 
would read 

dH 1 



(2tt) 4 (£ 2 -m 2 t + ie)(£ 2 - m 2 g + ie)(P - m| + ie) 



(18) 



A detailed calculation of this integral is given in the appendix. 
Using the above one finds for S(f)jk the contribution: 



(1) _ ^ g-mt 

where 



2^, 2 2 



5 nk = ^ySE^MJ^i^, - K b U k2 D b Jf(m 2 t ,ml,ml) (19) 
iD7r i=i i=i 



and 



1 z x y 

fix, y, z) = - r x (zxln- + xyln- + yzln-), (20) 

[x - y){x - z){z - y) x y z 

1 y 

f(x,y,y) = - r- x (xln- + x - y) (21) 

(y — x) 1 x 

Similarly one finds for the correction Aipjk from the same loop the following con- 
tribution 

A ^ = ^EEWW - ^U; 2 DtJf(mlml,ml) (22) 
i07r j=i i=i 



Next for the loop Fig. lii(a) we find 



A^ ( ? = (23) 
For the loop of Fig. li(b) we find 

tyfi = ^EE f «HW " K t V j2 D t2l )f(mlml,mD 
i07r i=i i=i 

2™ 2 2 



= Kft ifi^ S E F u U h D tu( V ki D *t u ~ KtV: 2 D* t Jf(mlmlml) (24) 
i07r i=i j=i 

where Fjj is given by 

F h = - ° Mz ((1 - ^ sin 2 6 w )D* tll D m + ^ sin 2 6 w D* m D t2i ) cos/5 
V 2 costly z o d 

» m '" :A - llA2 , (25) 



V^mvi/ sin (3 
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For the loop of Fig. lii(b) we find 



= 



2^,2 



A ^ = -^E^AW - KtV: 2 Dl)f(mlmlrni) (26) 



16vr 2 



i=i 



For loop of Fig. lii(c) we find 



4 4 

45 = 2 sEE #«4 sin /H- cos P 

i=l i=l 

m X? m X? „ 2 2 2 x 

16?r 2 /Ko,m x? ,m g+ ) 

Alrff = (27) 



where e and e are given by 



€ji - -gX^l - ~^ X 2i V i2 ~ ^ tan ^Xij-V^j 



e ;, = -^*.C/ a + ^^-^2 + ^= tanOwX^Ua (28) 



The parameters are defined as: 

<&■ = - tan^X^.)] (29) 

The matrix elements X are defined as 

X T M x oX = diag(m x o, m x o,m x o,m x o) (30) 
For loop of Fig. li(c) we find 



gm w cos (3 . 2 2 * 

P — [1 + 2 sin p — cos 2p tan uw\£^ 

2 v2 i=1 

4 sin ^ cos /3 

m X- 2 2 2 \ 

— \j [m ,m H+ ,m H+ ) 

(6 ) _ gm w cos/3 



$ = [1 + 2 sin 2 (3 - cos 2/3 tan 2 9 W \ £ 

2 v2 i= i 



tik cos p^* sin /3 



(31) 
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For loop of Fig. li(d) we find 



= ^ Q m /| CQS f E E E((^i - i^ms smP)(3Y n + iY l3 sin/5) 
«V^cos(% l=lm=li=1 

-2(Y m2 - %Y m3 cos/5)(y 2 + iY l3 cos/3) - AY m2 (Y n - iY l3 sin (3) tan/5) 

(Qki(Y n + sin/5) + S kl (Y l2 + «y I3 cos/5)) 
(Qi^mi + iY m3 sin/5) + S^(y m2 + iY m3 cos/5)) 

W% = 9 \ m r° S n E E E(0^i " sin/5)(3y u + iY l3 sin/5) 

»V^COS(% Hm=l i=l 

-2(r m2 - iy m3 cos/5)(y 2 + iY l3 cos/5) - 4Y m2 (y 1 - iY l3 sin/5) tan/5) 

(g* fc (Fa - iY l3 sin/5) + ^(y 2 - iY l3 cos/5)) 

(Q^(y m i - fy^a sin/5) + s;,(y m2 - *y m3 cos/5)) 

771 + 

^/K+.m^.m^,) (32) 

where Qji = an d Sji = -j^ipij, and the matrix elements Y are defined as 

YM 2 Higgs Y T = diag{m 2 H o,m 2 H o,m 2 H o). 
For loop of Fig. lii(d) we find 

$<f> { $ = -g 2 EEE4 

m=l i=l /=1 

(Q/i(y m i + ^y m 3 sin/5) + S^(y m2 + %Y m3 cos/5))(Q fci (y m i + «y m3 sin/5) 

m x +m x + 

+S ki {Y m2 + %Y m3 cos /5) ) * 6 ' / (m x + , % o , m x + ) 

AtfjJ = (33) 

For loop of Fig. lii(e) we find 



3 

2 2 



,(9)_ V 



*4% - ^EEW^,«0 P4) 



J=l i=l 

The parameters L' and R' are defined by 

1 



L l3 = -VaV^ - -y i2 y/ 2 + S tj sin 2 W 



4 = -Ut-yUjy - l -U* 2 U 32 + Sij sin 2 ^ (35) 
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For loop of Fig. li(e) we find 

(io) _ V2g 3 m z cosp ^ , , "< x , 2 2 , 



A4T = ^ (36) 



J COS^tV i^i 

For loop of Fig. lii(f) we find 



5& = 



4 4 



A^ X) = -4^ 3 E E ™y (37) 



167r 2 

=1 z=l xu/l 



where L, i? and Q" are defined as 

1 



2 

For loop of Fig. li(f) we find 



Rij = -^XziUj2 + X2iUj\ 

9Q" = \(X* 3i (gX* 2j - g'Xl ) + j)) (38) 



,,(12) _ *g 3 m w cos (3 * rn x o 

jk ~ ^2 ^ ij 16t^ ^ x * w+ ' w+ > 

A / (i2)_ 4g 3 m w cos (3 * m x o 

^Wjk ^ ^ttijUkj^J{m x o,m w+ ,m w +) {6\)) 

For loop of Fig. lii(g) we find 

m 3) = o 

A^ 3) = - g ^h T K T U; 2 V k \^f(m 2 T , m 2 T , r<) (40) 

where 

«r = 7= T (41) 

V^m^ cos p 

The loop corrections for and A^ fc are given by 

n=l 

A^ jk = £ A^ (42) 
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2.2 Loop analysis of A(j)jk and 5ipjk 

We do the same analysis of Fig. 2 as for Fig. I. We write down here the final 
results for both corrections from the thirteen loops together. The corrections are 
written in the same order of the loops in Fig. 2. 

= «t^EEW,W - K b U* 2 D* b2i )f(m 2 t ,ml,ml) 
iU7r i=i i=i 

2 2 2 

-Ktht^E v :2D bli (U;iD* bli - K b U; 2 DtJf(mlmlml) 
+K b j^j:j:EuU; 2 D tl[ (V k \D* tu - K t V: 2 Dl)f(m 2 b ,ml,ml) 



+0 
+0 

+ gm W Bm/3 ^ + 2 ^ ^ + ^ ^ * 

2 v2 i=1 

€ik cos /9ej* sin /3 

m X° 2 2 2 \ 

Y^/K°, %+-%+) 

+/ Q m ^ CQS f EE E(tan/3(F /2 -^ 3 cos/?)(3F m2 + ^ m3 cos/3) 

-4F n (F m2 - ir m3 cos /3) - 2 tan /3(y rol - iY m3 sin /3) (F a + iY l3 sin /3) ) 

(Q* k (Y n - iY l3 sin (3) + S* k (Y l2 - iY l3 cos/3)) 
{Q* 3l {Y ml - iY m3 sin/5) + S* 3l (Y m2 - iY m3 cos/5)) 

771 + 

Xa r/ 2 2 2 \ 

+0 



4g 2 2 2 m x +m X + 

cos 2 6 w ^^ ]l lk 16tt 2 jk z °' xt' 

V2g 3 m z sm(3 2 m + 

2^ R ji L ikJ7-^J K m x + > m Z > ™Z ) 



cos 3 # w ^ Jl " s 16tt 2 "' v *r 

4 4 ^2 0777, o 

i=i /=i iC)7r 

4<rWsin/5 ^ m x o 

^= 2^ Y67T 2 ^ ^ X * ' W+ ' W+ ' 



+0 (43) 
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where G and E are given by 

Gi i = J MZ q ((~5 + I sin2 °w)D* bli D blj - \ sin 2 9 w D* b2i D b2j ) sin/3 

+ n L> bli L> b2j 

\J lm w cos p 

E v = F> M \ ~ \ sin2 M^iiAy + I sin 2 ewD^Daj) sm/3 
V 2 cos 6V id o 



( D m D tij + a*2,a 2 ,) - » : - D* 2i D t2 j 



\ , lm w sm.p y/zmwsm.p 
and 5"' is given by 

The corrections Siftjk are given by 
2 2 



iD7r i=i ;=i 

+0 

2 2 2 

+ K SEE^M ll (W i - K t V j2 D t2l )f(mlmi,ml) 

i07r i=i ;=i 

+0 

4 4 

_2 f E E sin P e *j cos # 

i=l 1=1 

m ^ m xU- ( 2 2 2 , 



16tt 2 jx 

+ 2 cos 2 /3 + cos 2/3 tan 2 W ] V 
2V2 i=1 

4 sin ^ cos /3 



%f( m lo, m H+, m H+) 



167T 2 

777- 7 COS /3 ^ ^ ^ 

E E E( tan /^2 - ^3 cos/3)(3F m2 + iY m3 cos/3) 

8V2cos6V , =lm=u=1 

-4YJ 1 (y m2 - iVma cos/3) - 2tan/3(F m i - «y m3 sin/3)(Y"n + iY l3 sin/3)) 

+ sin/3) + S^Q^ + ^3 cos/3)) 
(Qii(V m i + ?F m 3 sin/3) + Sij(Y m2 + iY m3 cos/3)) 

m xt ci 2 2 2 \ 
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3 2 2 
m=l i=l 1=1 

(Qij(Y m i + iY m3 sin/5) + Sij(Y m2 + iY m3 cos j3)) (Q k i(Y ml + iY m3 sin/3) 

m x +m x + 

+S ki (Y m2 + iY m3 cos/3)) * 1 f(m +,m H0 ,m +) 

1071" m *I 

+0 

y^g 3 m z sin^ ™ x + 

+0 

4gWsin^ A m x o 
^= 1. L ij H ^Y^I ( m x ° ' m w+ , "V+ ) 

+0 (46) 



where S' is given by 



Sj,. = - tan 6^)] (47) 



3 Neutral Higgs decays including loop effects 

We summarize now the result of the analysis. Thus C e ff of Eq.(5) may be written 
as follows 

C eff = H?4(a% + Tscffixt + H.c (48) 

where 

af k = ^((>n+^3sin/3)(0 J - fc +<^ 

(49) 

and where 

ajf = ;^((^i+^3sm/3)(0 jfe +^ 

(50) 

Next we discuss the implications of the above result for the decay of the neutral 
Higgs. 



([kl«f,l) 2 + (l<l) 2 )(M^o - mj+ - m* + ) - \{{\a%\f - (|<|) 2 )(2m x +m x +)])(51) 
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There are many channels for H® decays. The important channels for the decay 
of the neutral Higgs boson are bb, it, ss, cc, fr, xtxj an d X^x'j- There is another 
set of channels that neutral Higgs can also decay into: these are modes of decaying 
into the other fermions of the SM, squarks, sleptons, other Higgs bosons, W and 
Z boson pairs, one Higgs and a vector boson, 77 pairs and finally into the gluonic 
decay i.e, Hf — > gg. We neglect the lightest SM fermions for the smallness of their 
couplings. We choose the region in the parameter space where we can ignore the 
other channels which either are not allowed kinematically or suppressed by their 
couplings. Thus in this work, squarks and sleptons are too heavy to be relevant 
in neutral Higgs decay. The neutral Higgs decays into nonsupersymmetric final 
states that involve gauge bosons and/or other Higgs bosons are ignored as well. 
In the region of large tan/5, these decays typically contribute less than 1% of the 
total Higgs decay rate [25]. Thus we can neglect these final states. 

We calculate the radiative corrected partial decay widths of the important 
channels mentioned above. In the case of CP violating case under investigation 
we use for the radiatively corrected V of neutral Higgs into quarks and leptons the 
analysis of [16], for the radiatively corrected partial widths into charginos we use 
the current analysis, and for the radiatively corrected decay width into neutralino 
we use [26]. We define 



where the first term in the numerator is the decay width including the full loop 
corrections and the second term is the decay width evaluated at the tree level. 
Finally to quantify the size of the loop effects on the branching ratios of the 
neutral Higgs decay we define the following quantity 



where the first term in the numerator is the branching ratio including the full loop 
corrections and the second term is the branching ratio evaluated at the tree level. 
The analysis of this section is utilized in Sec. (4) where we give a numerical analysis 
of the size of the loop effects and discuss the effect of the loop corrections on decay 
widths and branching ratios. 




(52) 




(53) 
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4 NUMERICAL ANALYSIS 



In this section we discuss in a quantitative fashion the size of loop effects on the 
partial decay width and the branching ratios of the neutral Higgs bosons into 
charginos. The analysis of Sec. 2 is quite general and valid for the minimal su- 
persymmetric standard model. For the sake of numerical analysis we will limit 
the parameter space by working within the framework of the SUGRA model [14]. 
Specifically we will work within the framework of the the extended mSUGRA 
model including CP phases. We take as our parameter space at the grand unifica- 
tion scale to be the following: the universal scalar mass m , the universal gaugino 
mass m 1 / 2 , the universal trilinear coupling |A |, the ratio of the Higgs vacuum ex- 
pectation values tan/3 =< H 2 > / < Hi > where H 2 gives mass to the up quarks 
and Hi gives mass to the down quarks and the leptons. In addition, we take for CP 
phases the following: the phase 9^ of the Higgs mixing parameter //, the phase ctA 
of the trilinear coupling A and the phases = 1,2,3) of the SU(3)c, SU(2) L 
and U (l)y gaugino masses. In this analysis the electroweak symmetry is broken by 
radiative effects which allows one to determine the magnitude of \i by fixing M z . 
In the analysis we use one loop renormalization group (RGEs) equations for the 
evolution of the soft susy breaking parameters and for the parameter /i, and two 
loop RGEs for the gauge and Yukawa couplings. In the numerical analysis we com- 
pute the loop corrections and also analyze their dependence on the phases. The 
masses of particles involved in the analysis are ordered as follows: for charginos 
m x + < m x + and for the neutral Higgs (m Hl ,m H2 ,mH 3 ) — > (jnH, r nh, m A) in the 
limit of no CP mixing where m H is the heavy CP even Higgs, m h is the light CP 
even Higgs, and itla is the CP odd Higgs. 

We investigate the question of how large loop corrections are relative to the 
tree values. We first discuss the magnitude of the loop corrections of the partial 
decay width defined in Eq.(52). As we mentioned earlier the loop corrections to 
the partial decay width of the chargino channel have been investigated before in 
the CP conserving case [20, 21]. The correction in these analyses is of the order of 
~ 10% of the tree level value. Our analysis supports this conclusion. In Figs. (3) 
and (4) we give a plot of Ar^ 1 ' 1 (/ = 1,3) as a function of tan/3 for the specific set 
of inputs given in the captions of these figures. We notice that the partial decay 
width gets a change of 7 ~ 15% of its tree level value. We also notice that the 
CP violating phase 9^ can affect the magnitude of this change. This effect has not 
been addressed in the previous analyses as they are working in the CP conserving 
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scenario. To compare between our analysis and the previous ones we have to notice 
that these analyses are using the general SUSY parameter space where they put 
by hand all the parameters that control the analysis. In [20], the authors choose 
the SUSY parameter set SPSla of the Snowmass Points and Slopes as a reference 
point. They choose for the trilinear couplings the values of A t = —487 GeV, 
Ab = —766 GeV and A T = —250 GeV. The values of the other parameters are: 
M = 197.6 GeV, M' = 98 GeV, fi = 353.1 GeV, tan/3 = 10, m A o = 393.6 GeV, 
Ms, = 558.9 GeV, M fj = 540.5 GeV, Mf, = 538.5 GeV, M f = 197.9 GeV, 

Vl,2 ' (A, 2 ' J->1,2 ' M,2 

M Sl 2 = 137.8 GeV, M Qa = 512.2 GeV, = 432.8 GeV, = 536.5 GeV, 
M La = 196.4 GeV and M^ 3 = 134.8 GeV. In all the figures of [20], these values 
are used, if not specified otherwise. In our mSUGRA analysis the magnitude of all 
these parameters and others are fixed by the five input parameters m = 100 GeV, 
m.i/2 = 250 GeV, tan/3 = 10, A = —100 GeV and a positive sign of fi in the CP 
conserving scenario [27]. These parameters are different from those of our Figs. 
(3) and (4). By using these parameters and fixing some of them by hand when 
needed to match their values in the analysis of [20], we were able to have a fair 
agreement with their Figs. (2-9). As an example of this check we show in Table. 1 
a comparison of the two works. For the input of Fig. 2 of [20] with CP violating 
phases are set to zero we can see that partial decay widths in both works have 
the same behavior as functions of masses and their magnitudes are fairly close to 
each other. However it seems that our loop corrected values of the partial widths 
are different from those of Eberl et al. This could be understood since our loop 
analysis of the effective lagrangian includes only the vertex corrections beside the 
corrections in the Higgs potential. 



case 


-ptree 
eberl 


ptree 
our 


-ptoop 

eberl 


-ploop 

our 


2.a m Ao = 700 GeV 


0.95 GeV 


0.94 GeV 


0.85 GeV 


0.80 GeV 


2.a m Ao = 800 GeV 


1.18 GeV 


1.17 GeV 


1.0 GeV 


0.91 GeV 


2.b m Ho = 800 GeV 


0.7 GeV 


0.69 GeV 


0.63 GeV 


0.58 GeV 


2.b m Ho = 900 GeV 


0.8 GeV 


0.8 GeV 


0.73 GeV 


0.70 GeV 



Table 1: A comparison between the current analysis and Eberl et al [20] for bench- 
mark cases. 

In the work of Ref. [21] only 8 out of 26 diagrams of the present analysis 
are calculated and they correspond to the vertex corrections from Figs. (l,2ii(a)), 
(l,2ii(b)), (l,2i(b)) and (l,2i(a)). By considering these diagrams only in the com- 
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parison, our analysis is in fair agreement with their Figs (2-4) and Figs. (6,8) for 
their inputs. 

Now we turn to address the question of how much loop corrections can affect 
the branching ratios into charginos. The branching ratio of a decay mode is defined 
to be the ratio between the partial decay rate of this mode and the total decay 
rate. In the parameter space under investigation this total decay rate includes the 
rates of decays into charginos, heavy quarks, taus and neutralinos. In Figs. (5) 
and (6) we give a plot of ABrj' 1 ^ = 1, 3) defined by Eq.(53) as a function of tan/3 
for the specific set of inputs given in the captions of these figures. Fig. (5) is for 
the neutral Higgs Hi boson and Fig. (6) is for the neutral Higgs H 3 boson. In all 
regions of the parameter space investigated in this work, the decay of the lightest 
Higgs boson H 2 into charginos is forbidden kinematically, since we have in these 
regions the fact that 2m^- > m# 2 . The analysis of Figs. (5) and (6) shows that 
the loop correction varies strongly with tan/5 with the correction changing sign 
for the case of H 3 decay. Further, the analysis shows that the loop correction can 
be as large as about —40% of the tree contribution for both Hi and H 3 cases. We 
also notice that the behavior of ABrj' 1 ^ = 1,3) as a function of tan/3 changes 
considerably by changing the phase of ji. So for some values of this phase we 
find that this parameter increases as tan/3 increases and for other values of 9^ 
we see that it decreases as tan (3 increases. As shown in the previous figures, the 
parameter tan (5 is playing a strong role. This parameter is important at the tree 

level through the diagonalizing mass matrices of the chargino and neutral Higgs 

IP s 

and their spectrum. At the loop level it has extra effect explicitly in a^ k ' and 
implicitly through the radiatively corrected matrix elements Y\ m and through the 
corrections 5(f) jk , A(f> jk , 5ipj k , Aif)j k . The values of the branching ratios themselves 
at tree and one loop levels are shown in Table. 2. 





Br°(Hi) 


Br lo °P(Hi) 


Br°(H 3 ) 


Br toop (H 3 ) 


0.5 


6% 


4.7% 


18.2% 


13.8% 


1.0 


8.4% 


6.9% 


21.3% 


18.1% 


1.5 


9.2% 


7.9% 


23.4% 


22.2% 



Table 2: Values of branching ratios at tree and one-loop levels of neutral Higgs 
into the channel XiXi at tan (3 = 24 for the input of Figs. (5) and (6) 

We notice that their magnitudes are not negligible for the region of the param- 
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eter space investigated. These non negligible branching ratios for the decay of the 
neutral Higgs into charginos suggest that these decay modes could be measurable 
at the soon-to-operate LHC. However, one should also consider the production 
rates for Hi and H$ bosons to assess whether the change in branching ratios could 
be detectable at colliders. This analysis goes beyond the scope of the current work. 
We also notice that the phase of the parameter \x affects the tree level branching 
ratios as well. This comes mainly from the structure of the chargino matrix. The 
more important channels in the region of the parameter space investigated are the 
decay into bottom and top quarks. They have the highest values of branching 
ratios. The radiative corrections of these channels are also more than those of 
the charginos and neutralinos. These channels were studied before [1, 15, 16] as 
mentioned above. However a 20% of branching ratio for the case of neutral Higgs 
as shown in the above table is not very small and could justify carrying out the 
current analysis. 

In Figs. (7) and (8) we give a plot of ABr]' l (l = 1,3) as a function of \A \ 
for the specific set of inputs given in the caption of these figures. The analysis of 
these figures shows that the loop corrections are substantial and reaches the value 
of —38% of the tree contribution for the case of Hi decay and the value of —43% 
for the case of H 3 decay. 

Next we investigate the effects of CP violating phases on the loop corrections 
of the neutral Higgs decays into charginos. In Figs. (9) and (10) we give a plot 
of ABr]' 1 ^ = 1,3) as a function of 9^ for the specific set of inputs given in the 
caption of these figures. The analysis of the figures shows that the loop correction 
has a sharp dependence on 9^. Further, the correction is changing sign as 9^ varies 
from to 7r for two cases of H 3 decay. Thus 9^ affects not only the magnitude of 
ABr]' 1 but also its sign depending on the value of 9^. 

In Figs. (11) and (12) we give a plot of ABr] ,l {l = 1,3) as a function of aA 
for the specific set of inputs given in the caption of these figures. Here also we find 
a very substantial dependence of ABrj' 1 on cha - This dependence is very large in 
the case of H 3 decay and it exceeds —40% of the tree contribution. 

In Figs. (13) and (14) we give a plot of ABrj' 1 ^ = 1, 3) as a function of £ 2 for 
the specific set of inputs given in the caption of these figures. Here we find a small 
effect of this phase on the loop corrections. 
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5 CONCLUSION 



In this paper we have carried out an analysis of the supersymmetric loop correc- 



breaking of electroweak symmetry one is left with three neutral Higgs bosons 
which in the absence of CP phases consist of two CP even Higgs bosons and one 
CP odd Higgs boson. In the absence of loop corrections, the lightest Higgs bo- 
son mass satisfies the inequality rrih < Mz and by including these corrections the 
lightest Higgs mass can be lifted above M z . With the inclusion of CP phases 
the Higgs boson mass eigenstates are no longer CP even and CP odd states when 
loop corrections to the Higgs boson mass matrix are included. Further, inclusion 
of loop corrections to the couplings of charginos and neutral Higgs is in general 
dependent on CP phases. Thus the decays of neutral Higgs into charginos can be 
sensitive to the loop corrections and to the CP violating phases. The effect of the 
supersymmetric loop corrections is found to to be in the range of 7 ~ 15% for the 
partial decay width. For the branching ratios it is found to be be rather large, as 
much as 40% in some regions of the parameter space. The effect of CP phases on 
the modifications of the partial decay width and the branching ratio is found to 
be substantial in some regions of the MSSM parameter space. 
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tions to XjXk H i 



couplings within MSSM. In supersymmetry after spontaneous 



6 APPENDIX 



The integral of import to this work is 



d 4 k 1 



(54) 



(2tt) 4 (Jfe 2 -m\ + ie)(k 2 -m 2 2 + ie)(k 2 - m\ + ie) 



It could be written in the form 




(55) 



where 



1 _ 1 1 1 

D ab c 
a = k 2 — ml + ie 
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b = k 2 — m\ + ie 

c = k 2 — m\ + ie (56) 

Using Feynman parametrization, could be written as 
1 f 1 f 1 - x 1 

— = 2 / dx dz-. - r -. --pr (57) 

D Jo Jo [a + (b - a)x + (c - a)z} 3 v ' 

The denominator in the above integral could be written in the form k 2 + M 2 + ie 
where M 2 = (ml — m 2 )x + (ml — ml)z — ml. Thus the integral J can take the 
form 

r d A k r 1 r l ~ x 1 

J = / - — -:2 / dx / dzj— - (58) 

J (2vr) 4 Jo Jo [k 2 + M 2 + ie] 3 V ; 

Now integrating over k and using the standard integral, for n > 3 

I = i7r2 r(n-2) j_ 

(2tt) 4 (P + A + ie)" T(n) A n ~ 2 w; 

one can find that the integral J has the form 

i 



J 



I dx dz — (60) 

(47r) 2 Jo Jo a + (3z 



where a = (ml — m 2 ):r — ml and /3 = m 2 — m 2 . Integrating over z one can get for 
the integral J the form of 

J = 



i 1 f 1 

, - xo — ^ 5-/ dxln(5ix — m\) — ln(<5 2 x — m?) (61) 

(47r) 2 mf - m§ Jo 

where <5i = m 2 — m 2 and 82 = ml — m\. Finally we integrate over x to get for J 
the form of 

J = j^yf(mlm 2 2 ,m 2 3 ) (62) 



where 



fi 2 2 2\ _ *■ J - 

j {m 1 ,m 2 ,m 3 ) — „- 2 2 2 2 2 
- WI3 WI3 — 7712 mf — wij 

222 

x \mlm\ ln(^|) + mlm 2 , ln(^f ) + m?m? ln(— ?)1 (63) 
L 2 d v m§ y 6 1 v m 2/ 1 2 m\ n v 7 

This is the famous form factor that appears in the analysis of the radiative correc- 
tions for the quark and lepton masses [28] , the decay rates of neutral and charged 
Higgs into quarks and leptons [1, 15, 16, 18] and in the b — > 57 process [17]. In the 
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2 

latter process, the authors are using different form factor This form 

factor could be easily converted to our f (m\, m\, m 2 ,) through the simple relation 



2 2 

r mi mi 



mlf{mlmlml)=HQ^) (64) 

For the case where two of the masses are equal, m 2 = m 3 , one can repeat the 
same analysis with b = c = k 2 — m\ + ie and a = k 2 — m\ + ie. By doing so one 
can get for the form factor J 

% 1 r 2 i / ^^"^ \ 9 2n / ^ \ 

(47r) 2 (m% - m{) 2 m{ 
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Figure 1: Set of diagrams contributing to radiative corrections 5<pjk and Aipj^. (i): 
(a) si = b*, s 2 = bf, f = t; (b) s x = t h s 2 = t h f = b; (c) s x = H+, s 2 = H+, 

f = xl (d) si = fl?, s 2 = / = xt; (e) «i = z°, ^ 2 = z°,f = xt; (f) 

= W+, t; 2 = W+, f = xl (ii): (a) fx = t, f 2 = t, s = b*; (b) /i = 6, / 2 = b, 
s = k (c)/i = X°, /a = X?, s = # + ; (d) fx = xt, h = xt,s = H° m , (e) f\ = xt, 
f 2 = xt, v = Z°; (f) fx = x°i, / 2 =Xi,v = W- (g) /i = r+, / 2 = r+, s = v T . 
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Figure 2: Set of diagrams contributing to radiative corrections A0 jfc and Sifrjk- (i): 
(a) si = b*, s 2 = bf, f = t; (b) s x = t h s 2 = t h f = b; (c) s x = H+, s 2 = H+, 
f = xl (d) si = Hi s 2 = H° m , f = X f; (e) v, = Z°, v 2 = Z°, f = xt; (f) 
v x = W+, v 2 = W + ,f = xl (ii): (a) fi = t, f 2 = t, s = b*; (b) h = b,f 2 = b, 
s = k (c)/i = X°, /a = X?, s = H+- (d) /i = X + f 2 = xt,s = H° m , (e) A = X + 
/ 2 = x + i, = (f) /! = xl / 2 = X?, « = (g) £ = r+, / 2 = r+, s = v T . 
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Figure 3: tan/3 dependence of Ar\ — > Xi Xi ■ The curves in ascending order 
correspond to # M = 0.2, 0.4, 0.6 (rad). The input is m = 350 GeV, mi/ 2 = 180 
GeV, 6 = 0.4 (rad), £ 2 = 0.5 (rad), £ 3 = 0.6 (rad), a Ao = 0.8 (rad) and \A \ = 250 
GeV. 



24 




Figure 4: tan f3 dependence of Ar 3 — * XiXi ■ The curves in ascending order 
correspond to 6^ = 0.2, 0.4, 0.6 (rad). The input is m = 350 GeV, mi/ 2 = 180 
GeV, £i = 0.4 (rad), £ 2 = 0.5 (rad), £ 3 = 0.6 (rad), a Ao = 0.8 (rad) and \A \ = 250 
GeV. 



25 




Figure 5: tan (3 dependence of ABri — > . The curves in ascending order 

at tan/? = 40 correspond to 9^ = 0.5, 0.1, 1.0, 1.5 and 2.0 (rad). The input is 
m = 500 GeV, m 1/2 = 150 GeV, & = 0.4 (rad), & = 0.5 (rad), £ 3 = 0.6 (rad), 
a Ao = 0.3 (rad) and \A \ = 250 GeV. 
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Figure 6: tan j3 dependence of ABr 3 — > . The curves in ascending order 

at tan/? = 40 correspond to 9 ^ = 0.5, 0.1, 1.0, 1.5 and 2.0(rad). The input is 
m = 500 GeV, m 1/2 = 150 GeV, & = 0.4 (rad), £ 2 = 0.5 (rad), £ 3 = 0.6 (rad), 
a Ao = 0.3 (rad) and \A \ = 250 GeV. 



27 




CD 

o 



Figure 7: |A | dependence of ABri — > Xi~Xi • The curves in ascending order at 
\A \ = correspond to tan/5 = 40, 35, 30, 25 and 20. The input is m = 500 GeV, 
mi/j = 150 GeV, & = 0.4 (rad), £ 2 = 0.5 (rad), £ 3 = 0.6 (rad), M = 0.7 (rad) and 
tt^o =0.1 (rad). 
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Figure 8: \A \ dependence of ABr 3 — > XiXi ■ The curves in ascending order at 
\Aq\ = correspond to tan/5 = 40, 35, 30, 25 and 20. The input is m = 500 GeV, 
mi/j = 150 GeV, & = 0.4 (rad), £ 2 = 0.5 (rad), £ 3 = 0.6 (rad), M = 0.7 (rad) and 
®a =0.1 (rad). 
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Figure 9: 9 ^ dependence of ABr 1 — > XiXi ■ The curves in ascending order at 
9^ = 2.0 (rad) correspond to \A \ = 100, 250, 500, 750 and 900 GeV. The input 
is tan/5 = 20.0, m = 500 GeV, m 1/2 = 150 GeV, & = 0.4 (rad), £ 2 = 0.5 (rad), 
£ 3 = 0.6 (rad) and cxa = 0.2 (rad). 
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Figure 10: 9 ^ dependence of ABr 3 — > XiXi ■ The curves in ascending order at 
9^ = n (rad) correspond to \A \ = 100, 250, 500, 750 and 900 GeV. The input is 
tan/5 = 20.0, m = 500 GeV, m 1/2 = 150 GeV, & = 0.4 (rad), & = 0.5 (rad), 
£ 3 = 0.6 (rad) and cxa = 0.2 (rad). 
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Figure 11: «o dependence of ABr\ — > XiXi ■ The curves in ascending order at 
a Ao = 2.2 (rad) correspond to \A \ = 500, 450, 400, 100 and 200 GeV. The input 
is tan/5 = 20.0, m = 500 GeV, m 1/2 = 150 GeV, = 0.4 (rad), £ 2 = 0.5 (rad), 
£ 3 = 0.6 (rad) and 9 ^ = 0.1 (rad). 
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Figure 12: a dependence of A_Br 3 — > Xi Xi • The curves in ascending order at 
a Ao = 2.2 (rad) correspond to \A \ = 500, 450, 400, 100 and 200 GeV. The input 
is tan/5 = 20.0, m = 500 GeV, m 1/2 = 150 GeV, = 0.4 (rad), £ 2 = 0.5 (rad), 
£ 3 = 0.6 (rad) and 9 ^ = 0.1 (rad). 
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Figure 13: £ 2 dependence of ABri — > XiXi ■ The curves in ascending order at 
£ 2 = 0.75 (rad) correspond to \A \ = 50, 100, 150, 200 and 250 GeV. The input is 
tan (3 = 20.0, m = 500 GeV, m 1/2 = 150 GeV, & = 0.4 (rad), £ 3 = 0.6 (rad) and 
6 ^ = 0.2 (rad) and cxa = 0.3 (rad). 
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Figure 14: £ 2 dependence of A5r 3 — > Xi Xi • The curves in ascending order at 
^ = 0.75 (rad) correspond to |A | = 50, 100, 150, 200 and 250 GeV. The input is 
tan (3 = 20.0, m = 500 GeV, m 1/2 = 150 GeV, & = 0.4 (rad), £ 3 = 0.6 (rad) and 
6> M = 0.2 (rad) and au = 0.3 (rad). 
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